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BESCHREIBUNG 

Computertomographie-Verfahren mit redundanten Messwerten 

Die Erfindung betiiffl ein Coir^utertomograpWe-Verfehren, bei dem ausschlieBlich 
Messwerte zurRekonstaiktion eines CT-Bildes verwendet werden, die in einem 
5 Rekonstruktionsfenster liegen, insbesondere betrifft die Erfindung exakte Computer- 
tomographie-Verfohren. Die Erfindung bezieht sich auBerdem auf einen Computer- 
tomographen zur Durchfuhrung dieses Verfahrens sowie auf ein Computerprogramm 
zur Steuerung des Computertomographen. 

10 In bekannten Verfahren der eingangs genannten Art bewegt sich eine Strahlenquelle 
relativ zu einem Untersuchungsbereich, in dem sich ein zu untersuchendes Objekt 
befindet, auf einer helix- oder kreisformigen Trajektorie. Dabei werden Messwerte mit 
einer Detektoreinheit akquiriert, aus denen sich ein CT-Bild, bspw. eine Absorptions- 
verteilimg im Untersuchungsbereich, rekonstruieten lasst Viele Rekonstruktionsver- 

15 fahren, insbesondere die exakten Rekonstruktionsverfahren, konnen nur einen Teil der 
Messwerte verwenden. Der Teil der Messwerte, der zur Rekonstruktion verwendet 
werden kann, kann bspw. dutch einen Raumwinkelbereich, in dem sich die Strahlen- 
quelle bewegt, wahrend eine zu rekonstruierende Stelle des Untersuchungsbereichs 
beleuchtet wird, oder durch einen bestimmten Bereich auf der Detektoreinheit definiert 

20 sein. Der Bereich, in dem sich Messwerte befinden, die zur Rekonstruktion verwendet 
werden, wird Rekonstruktionsfenster genannt Messwerte auBerhalb des Rekonstruk- 
tionsfensters, werden von diesen Verfahren nicht beriicksichtigt und als redundante 
Messwerte bezeichnet. Ein derartiges Rekonstruktions verfahren ist bspw. von Katsevich 
in 5r Analysis of an Exact Inversion Algorithm for Spiral Cone-Beam CT" Physics 

25 Medicine and Biology, vol. 47, pp. 2583-2597 (El) beschrieben worden. Dort bewegt 
sich die Strahlenquelle relativ zum Untersuchungsbereich auf einer helixformigen 
Trajektorie, und es werden nur Messwerte verwendet, die in dem sogenannten PI- 
Fenster, das unter naher erlautert wird, liegen. Nachteilig ist bei diesen Verfahren, dass 
sie ein schlechtes Signal-Rausch-Verhaltnis auftveisen und bei Objektbewegungen zu 

30 storenden Bewegungsartefiakten fuhren. 
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Es ist daher Aufgabe der vorliegenden Erfindung, ein Computertomographie-Verfahren 
der eingangs genannten Art anzugeben, bei dem das Signal-Rausch-Verhaltnis erhdht 
wird und die Bewegungsartefakte reduziert werden. 

5 

Diese Aufgabe wird erfindungsgemaB gelost durch ein Computertomographie-Ver- 
fahren mit den Schritten: 

a) Erzeugen eines kegeliormigen, einen Untersuchungsbereich und ein darin 
befindliches Objekt durchsetzenden Strahlenbundels mit einer Strahlenquelle, 
10 b) Erzeugen einer Relativbewegung zwischen der Strahlenquelle einerseits und 
dem Untersuchungsbereich andererseits, die mindestens eine 
Rotationsbewegung urn eine Rotationsachse umfasst und insbesondere die Form 
einer Helix oder eines Kreises hat, 

c) Akquirieren von Messwerten, die von der Intensitat in dem Strahlenbiindel 
15 jenseits des Untersuchungsbereichs abhangen, mit einer Detektoreinheit 

wahrend der Relativbewegung, 

d) Ermitteln von jeweils einem Komplementarmesswert fur zumindest einige der in 
Schritt c) akquirierten und in einem Rekonstruktionsfenster liegenden 
Messwerte mit Hilfe von redundanten Messwerten, wobei die dem jeweiligen 

20 Messwert und dem dazugehorigen Komplementarmesswert zugeordneten 

Strahlen entgegengesetzt zueinander orientiert sind, 

e) Ersetzen jedes Messwertes, fur den in Schritt d) ein Komplementarmesswert 
ermittelt worden ist, durch eine Summe aus gewichtetem Messwert und 
gewichtetem Komplementarmesswert, 

25 f) Rekonstruieren eines CT-Bildes des Untersuchungsbereichs aus den im 
Rekonstruktionsfenster liegenden Messwerten. 

Im Gegensatz zu bekaimten Computertomographie-Verfahren der eingangs genannten 
Art werden mit Hilfe von redundanten Messwerten Komplementarmesswerte ermittelt, 
30 die addiert auf die jeweils dazugehorigen Messwerte diesen ersetzen, so dass zu der 
folgenden Rekonstruktion auch redundante Messwerte beitragen. Die Erfindung beruht 
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nun auf der Erkenntnis, dass bei Computertomographie-Verfahren, die bei der Rekon- 
struktion nicht auf einen Teil der Messwerte beschrankt sind, Bewegungsartefakte 
reduziert sind und das Signal-Rausch-Verhaltnis erhoht ist Daher fuhrt die Beruck- 
sichtigung der redundanten Messwerte neben einer Erhohung des Signal-Rausch- 
5 Verhaltnisses zu einer Reduzierung der Bewegungsartefakte. 

Im Anspruch 2 wird der Komplementarmesswert, wenn er nicht in Schritt c) akquiriert 
worden ist, mit Efilfe von Johns-Gleichung ermittelt Dies fuhrt zu CT-Bildern 
besonders guter Qualitat. 

10 

Johns-Gleichung ist eine bekannte partielle Differentialgleichung, die es unter anderem 
ennoglicht, zu akquirierten Messwerten weitere Messwerte zu ermitteln, die mit den 
akquirierten Messwerten konsistent sind. 

15 Der Begriff Consistent" lasst sich wie folgt erMaxen. Durch den im Schritt c) akquirier- 
ten Satz an Messwerten ist das Objekt im Untersuchungsbereich nicht eindeutig fest- 
gelegt, da die Messwerte verrauscht sind oder bspw. auch zu wenig Messwerte akqui- 
riert worden sein konnen. Zu dem Satz an Messwerten gibt es also eine Menge an 
moglichen Objekten im Untersuchimgsbereich, die zu den Messwerten passen. Ein 

20 Komplementarmesswert ist konsistent zu den akquirierten Messwerten, wenn er durch 
Durchstrahlen eines der moglichen Objekte erzeugt werden kann. Eine Konsistenz- 
bedingung, mit der ein derartiger Komplementarmesswert ermittelt werden kann, ist 
bspw. die bekannte Johns-Gleichung. 

25 Im Anspruch 3 wird bei der Addition der gewichteten Komplementarmesswerte und der 
jeweils dazugehorigen gewichteten Messwerte ein akquirierter Komplementarmesswert 
starker gewichtet, wenn er im Schritt c) akquiriert und nicht bspw. durch Johns- 
Gleichung ermittelt worden ist Die starkere Beriicksichtigung wirklich gemessener 
Werte fuhrt zu einer verbesserten Badqualitat 
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Die Anspriiche 4 und 5 beschreiben ein bevorzugtes exaktes Computertomographie- 
Verfahren, dass fur helixformige Trajektorien zu CT-Bildem guter Qualitat fiihrt. 

Ein Computertomograph zur Durchfuhrung des erfindungsgemaBen Verfahrens ist in 
5 Anspruch 6 beschrieben. Anspruch 7 definiert ein Computerprogramm zur Steuerung 
des Computertomographen nach Anspruch 6. 

Die Erfindung wird im Folgenden anhand der Zeichnungen erlautert. Es zeigen: 

10 Fig. 1 schematisch einen Computertomographen gemaB einer besonderen Ausfiihrungs- 
form, mit dem das erfindungsgemafle Verfahren gemaB einer besonderen Ausfuhrungs- 
form ausfuhrbar ist, 

Fig. 2 ein Ablaufdiagramm des erfindungsgemaBen Verfehrens, 

15 

Fig. 3 eine PI-Gerade und ein Pl-Intervall fur einen Punkt im Untersuchungsbereich, 

Fig. 4 die PI-Gerade und das Pl-Ihtervall fur einen Punkt im Untersuchungsbereich 
projiziert in eine Ebene senkrecht zur Rotationsachse, 

20 

Fig. 5 schematisch einen planaren Detektor mit einem PI-Fenster, 

Fig. 6 eine schematische Schnittansicht einer helixfbnnigen Trajektorie und eines 
Untersuchungsbereiches mit einem direkten Strahl und einem Komplementarstrahl, 

25 

Fig. 7 ein Ablaufdiagramm fur eine exakte Rekonstruktion eines CT-Bildes 5 

Fig. 8 eine schematische Perspektivansicht paralleler Strahlen mit verschiedenen 
Strahlenpositionen, 

30 

Fig. 9 eine schematische Perspektivansicht der helixformigen Trajektorie mit einer k- 
Ebene und einer K-Linie, 
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Fig. 10 eine schematische Perspektivansicht der helixlSrmigen Trajektorie mit durch 
Rebinning gebildeten Strahlenfachern in paiallelen Ebenen, 

Der in Fig. 1 dargestellte Computertomograph umfesst eine Gantry 1 , die urn eine 
5 parallel zur z-Richtung des in Fig. 1 dargestellten Koordinatensystems verlaufende 
Rotationsachse 14 rotieren kann. Dazu wird die Gantry 1 von einem Motor 2 mit einer 
vorzugsweise konstanten, aber einstellbaren Winkelgeschwindigkeit angetrieben. An 
der Gantry 1 ist eine Strahlenquelle S befestigt, beispielsweise ein Rontgenstrahler. 
Dieser ist mit einer Kollimatoranordnung 3 versehen, die aus der von der Strahlenquelle 
10 S erzeugten Strahlung ein kegeliSrmiges Strahlenbiindel 4 ausblendet, d.h. ein Strahlen- 
biindel, das sowohl in z-Richtung als auch in einer dazu senkrechten Richtung (d.h. in 
einer zur Rotationsachse senkrechten Ebene) eine von Null verschiedene, endliche 
Ausdehnung hat 

15 Das Strahlenbiindel 4 durchdringt einen zylinderfonnigen Untersuchungsbereich 1 3, in 
dem sich ein Objekt, z. B. ein Patient auf einem Patientenlagerungstisch (beides nicht 
dargestellt) oder aber auch ein technisches Objekt befinden kann. Nach dem Durch- 
setzen des Untersuchungsbereichs 13 trifft das Strahlenbiindel 4 auf eine an der Gantry 
1 befestigte Detektoieinheit 16 mit einer Detektorflache, die eine Vielzahl von Detek- 

20 torelementen umfasst, die in dieser Ausfuhrungsform in Zeilen und Spalten matrix- 
formig angeordnet sind. Die Detektorspalten verlaufen parallel zur Rotationsachse 14, 
Die Detektorzeilen befinden sich in zur Rotationsachse senkrechten Ebenen, in dieser 
Ausfuhrungsform auf einem Kreisbogen urn die Strahlenquelle S (fokus-zentrierter 
Detektor). In anderen Ausfuhrungsformen konnen sie aber auch anders geformt sein, z. 

25 B. einen Kreisbogen urn die Rotationsachse 14 beschreiben oder geradlinig sein. Jedes 
von dem Strahlenbiindel 4 getroffene Detektorelement liefert in jeder Position der 
Strahlenquelle einen Messwert Sir einen Strahl aus dem Strahlenbiindel 4. 
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Der mit Omax bezeichnete Of&ungswinkel des Strahlenbundels 4 bestimmt den Durch- 
messer des Objektzylinders, innerhalb dessen sich das zu untersuchende Objekt bei der 
Akquisition der Messwerte befindet Dabei ist der Offiiungswinkel als der Winkel 
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definiert, den ein Strahl, der in einer zur Rotations achse 14 senkrechten Ebene am 
Rande des Strahlenbiindels 4 liegt, mit einer durch die Strahlenquelle S und die 
Rotationsachse 14 definierten Ebene einschlieBt. Der Untersuchungsbereich 13 bzw. das 
Objekt oder der Patientenlagerungstisch konnen mittels eines Motors 5 parallel zur 
5 Rotationsachse 14 bzw. zur z-Achse verschoben werden. Dazu aquivalent konnte aber 
auch die Gantry in diese Richtung verschoben werden. Wenn es sich urn ein technisches 
Objekt handelt und nicht urn einen Patienten, kann das Objekt bei einer Untersuchung 
gedreht werden, wahrend die Strahlenquelle S und die Detektoreinheit 16 still stehen. 

10 Wenn die Motoren 2 und 5 gleichzeitig laufen, beschreiben die Strahlenquelle S und die 
Detektoreinheit 16 eine helixformige Trajektorie relativ zum Untersuchungsbereich 13. 
Wenn hingegen der Motor 5 fur den Vorschub in Richtung der Rotationsachse 14 still- 
steht und der Motor 2 die Gantry rotieren lasst, ergibt sich eine kreisfiirmige Trajektorie 
fur die Strahlenquelle S und die Detektoreinheit 16 relativ zum Untersuchungsbereich 

15 13. In diesem Ausfuhrungsbeispiel wird nur die helixformige Trajektorie betrachtet. 
Das erfindungsgemaBe Verfahren ware aber auch bei einer kreisformigen Trajektorie 
anwendbar. 

Die von der Detektoreinheit 1 6 akquirierten Messwerte werden einem Rekonstruktions- 
20 rechner 10 zugefiihrt, der mit der Detektoreinheit 16 z. B. uber eine kontaktlos arbei- 
tende Dateniibertragting (nicht dargestellt) verbunden ist. Der Rekonstruktionsrechner 
10 rekonstruiert das CT-Bild und gibt es, beispielsweise auf einem Monitor 1 1, wieder. 
Die beiden Motoren 2 und 5, der Rekonstruktionsrechner 10, die Strahlenquelle S und 
der Transfer der Messwerte von der Detektoreinheit 16 zum Rekonstruktionsrechner 10 
25 werden von einer Steuereinheit 7 gesteuert 

Fig. 2 zeigt den Ablauf eines Mess- und Rekonstniktionsverfahrens, das mit dem 
Computertomographen nach Fig. 1 durchgefiihrt werden kann. 

30 Nach der Initialisierung im Schritt 101 rotiert die Gantry mit einer Winkelgeschwindig- 
keit, die in diesem Ausfuhrungsbeispiel konstant ist. Sie kann aber auch variieren, z.B. 
in Abhangigkeit von der Zeit oder von der Strahlenquellenposition. 
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Im Schritt 103 wird der Untersuchungsbereich bzw. das Objekt oder der Patienten- 
lagerungstisch parallel zur Rotationsachse verschoben und die Strahlung der Strahlen- 
quelle S wird eingeschaltet, so dass die Detektoreinheit 16 die Strahlung aus einer 
Vielzahl von Winkelpositionen erfassen kann. Die Strahlenquelle S bewegt sich dabei 
relativ zum Untersuchungsbereich auf einer helixformigen Trajektorie, die durch 



y(^) = 



'RcosX} 
hX 



(1) 



10 parametrisiert weiden kann. Hierbei ist 2%h der Tischvorschub pro Rotation, Z die 
Winkelposition der Strahlenquelle bezogen auf eine beliebige, aber feste Referenz- 
winkelposition und R der Radius der helixformigen Trajektorie 17. 

In diesem Ausfiihrungsbeispiel wird der Tischvorschub pro Rotation so gewahlt, dass 
15 von jedem Ort x im Untersuchungsbereich aus die Strahlenquelle S fiber einen 

Winkelbereich von zumindest 1 80° sichtbar ist. Die Aufhahme von Messwerten mit 
diesem Tischvorschub wird als PI-Akquisition bezeichnet. 

Die PI-Akquisition soli im Folgenden naher erlautert werden. Dazu sind in Fig. 3 die 
20 helixformige Trajektorie 17, auf der sich die Strahlenquelle relativ zur Stelle x im 
Untersuchungsbereich bewegt, und der Abschnitt / PI (x) der Helix 17, der durch eine 
PI-Gerade 37 eingeschlossen ist, dargestellt. Fig. 4 zeigt die Helix 17 aus Fig. 3 proji- 
ziert auf eine senkrecht zur Rotationsachse orientierten Ebene. Die PI-Gerade 37 ist die 
Linie, die die Helix an zwei Stellen xmd die Stelle x schneidet, wobei der durch die 
25 Gerade eingeschlossene Helixabschnitt 7 PI (x) einen Winkel kleiner 2n iiberstreicht. Bei 
der PI-Akquisition durchlaufen Strahlen, die von der Strahlenquelle S ausgehen, 
wahrend sie sich auf dem Helixabschnitt I n (x) befindet, die Stelle x. 



30 



Im Rahmen dieses Ausfiihrungsbeispiel wird unten ein Rekonstruktionsverfahren 
verwendet, dass ausschlieBlich Messwerte verwendet, deren zugeordnete Strahlen vom 
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Helixabschnitt 7 M (x) ausgehen, A h., es werden nur Messwerte verwendet, die im 
sogenannten Pl-Intervall liegen. Bekanntlich kann durch geometrische Betrachtungen 
gezeigt werden, dass Messwerte, deren Strahlen von einem bestimmten Helixabschnitt 
ausgehen, auf der Detektoreinheit in einem Rekonstruktionsfenster liegen. Bei der hier 
5 vorliegenden PI-Akquisition ist das Rekonstruktionsfenster das PI-Fenster. Das PI- 
Fenster 25 ist auf einem planaren Detektor 60 durch die PI-Linien 21, 23 begrenzt, die 
durch folgende Gleichungen definiert sind: 



15 Hierbei sind u und v Koordinaten auf dem planaren Detektor 60 gemaB dem Koordi- 
natensystem in Fig. 5. Dieses Koordinatensystem ist aus tJbersichtlichkeitsgrQnden 
unterhalb des planaren Detektors 60 dargestellt. Der Ursprung des Koordinatensystems 
liegt aber im Zentrum des Detektors. Die Einheitsvektoren dieses Koordinatensystems 
lauten l u = (-sink, cos k,0) und l v = (0,0,1) . 

20 

Der planare Detektor 60 ist ein fiktiver Detektor, der die Rotationsachse 14 enthalt und 
senkrecht zu dem Strahl orientiert ist, der ausgehend von der jeweiligen Strahlen- 
quellenposition senkrecht die Rotationsachse 14 triffi. Messwerte, die mit der realen 
Detektoreinheit 16 detektiert werden, konnen entlang der entsprechenden Strahlen auf 
25 den planaren Detektor 60 projiziert werden. 




(1) 



und 




(2) 



Computertomographie-Verfahren der eingangs genannten Art verwenden nur Mess- 
werte, die in einem Rekonstruktionsfenster liegen. Die ubrigen Messwerte werden nicht 
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benutzt, was zu einem schlechten Signal-Rausch-Verhaltnis und Bewegungsartefakten 
fuhrt. Ia den folgenden Schritten wird zu jedem Messwert g(X 0 ,u 09 v 0 ) 9 der im Rekon- 
struktionsfenster liegt, ein Komplementarmesswert g c (X 0 , u 0 , v 0 ) mit Hilfe von redun- 

danten Daten, also mit Messwerten, die auBerhalb des Rekonstruktionsfensters liegen, 
5 ennittel^ wodurch das Signal-Rausch-Verhaltnis veibessert und die Bewegungsarte- 
fakte reduziert werden. Dabei wird der dem Messwert g(X 0 ,u 0 ,v 0 ) zugeordnete Strahl 
im Folgenden direkter Strahl genannt 

In anderen Ausfuhrungsformen kann auch nur fur einen Teil der akquirierten Mess- 
10 werte, die im Rekonstruktionsfenster liegen, jeweils ein Komplementarmesswert 

bestimmt werden. Insbesondere ware es moglich, nur fiir diejenigen Messwerte einen 
Komplementarmesswert zu bestimmen, die in Randbereichen des Rekonstruktions- 
fensters liegen. 

1 5 Ein Komplementarmesswert ist durch den Verlauf des dazugehorigen Komplementar- 
strahles definiert. Ein Komplementarstrahl verlauft entgegengesetzt parallel zu dem 
dazugehorigen direkten Strahl. Zur Ermittlung des dem Komplementarstrahl zugeord- 
neten Komplementarmesswertes kann eine Konsistenzbedingung verwendet werden. 
Die Konsistenzbedingung, die in diesem Ausfuhrungsbei^piel zur Ermittlung des zum 

20 Komplementarstrahl gehorenden Komplementarmesswertes verwendet wird, ist die in 
„Partial Differential Equations", F. John, Applied Mathematical Sciences, Springer- 
Verlag, 1971 veroffentlichte allgemein mit Johns-Gleichung bezeichnete Gleichung, die 
in US 6,292,526 Bl fiir helixformige Trajektorien angepasst worden ist 

25 ErfindungsgemaB ist jedes Verfahren zur Bestimmung eines Komplementarmesswertes, 
insbesondere jede Konsistenzbedingung, anwendbar, das es ermogiicht, zu einem 
entgegengesetzt parallel zum direkten Strahl verlaufenden Komplementarstrahl einen 
Komplementarmesswert zumindest naherungsweise unter Berucksichtigung von 
redundanten Messwerten zu ermitteln. 

30 
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ZumEnnittelndesjeweiKgenKomplementarmessw £ c (^ 0 > tt o> v o) wirdzunachst 
im Schritt 105 aus den akquirierten Messwerten ein Naherungsmesswert g(X x ,u x ,v x ) 
bestimmt, dessen zugeordneter Naherungsstrahl 3 1 naherungsweise gleich verlauft wie 
der Komplementarstrahl (siehe Fig, 6). Dieser Naherungsstrahl geht von der 
5 Winkelposition X x aus und trifft auf den planaren Detektor 60 auf die Stelle (u x , v x ) . 

Die Winkelposition des Naherungsstrahles ist durch 



deflniert. Hierbei ist p der Facherwinkel des direkten Strahles, d.h., der Winkel, den der 
direkte Strahl mit einem Strahl einschlieBt, der den planaren Detektor 60 in Fig. 5 mittig 
schneidet Die Winkelposition X x ist die Winkelposition, an der der direkte Strahl, der 
von der Strahlenquellenposition y(X 0 ) ausgehend auf den planaren Detektor auf die 
1 5 Stelle (u 0 , v 0 ) trifft, erstmals die Helix 17 beruhrt, wenn er parallel zur Rotationsachse 
14 verschoben wird 

Fiir die i/-Koordinate des Naherungsstrahles auf dem planaren Detektors 60 gilt: 
20 u x = -u 0 . (4) 

Die v-Koordinate wird wie folgt ermittelt Durch X = X x und u = u x ist ein Strahlen- 
Scher definiert, der von der Winkelposition X } auf der Helix 17 ausgeht xmd auf alle 
Detektorelemente trifit, fur die u = u x gilt. Fiir jeden Strahl dieses Strahlenfachers wird 
25 der jeweilige Schnittpunkt mit dem direkten Strahl 33 ermittelt. Der Naherungsstrahl 3 1 
ist dann derjenige Strahl des Strahlenlachers, dessen Schnittpunkt am nachsten an der 
Rotationsachse 14 liegt. Dieser Naherungsstrahl trifift auf dem planaren Detektor auf die 
Stelle (u l9 v x ) mit 




X 0 +7i;-2p fur v 0 >0 
X 0 -n-2$ fur v 0 <0 



(3) 
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fif^R 2 



R 



A(2P— n) 



i <fif+ii*} 



2cosP 



R 



R 



*( 2 p +n )|_L__3g±£ 
A 2cosp i? 



-2v 0 cosp 



-2v 0 cosp 



fur v 0 >0 



fur v ft <0 



(5) 



Im Schritt 107 wind der Verlauf des Komplementarstrahles ermittelt. Der Komplemen- 
tarstrahl verlauft, wie oben bereits erwahnt, entgegengesetzt parallel zum direkten 
5 Strahl 33. Der direkte Strahl 33 trifft die Zylinderoberflache 35, die durch die helix- 
formige Trajektorie 17 definiert ist, an einer Stelle 



z = z 0 =hX 0 +v 0 



2J?cosp 



(6) 



10 Die Strahlenquelle S befand sich wahrend der Akquisition 
der Stelle 



an 



z = z 1 = hK l = 



/r(X, 0 +7t-2p) fur v 0 >0 
h(X 0 -n - 2p) fur v 0 < 0 



(7) 



15 Daher muss, urn einen Komplementarstrahl zu erhalten, die Strahlenqpiellenposition z x 
zur Position z 0 „verschoben" werden. Sie muss also urn 



Aq = z t -z 0 =■ 



_„ 2/2 cos P ~. _ 
-2B)-v n - F fur v 0 >0 



/i(+7t-2p) 



M-^t-2p)-v 0 



2R cosp 



(8) 



fur v„ <0 



20 verschoben werden. Die Position der Strahlenquelle fur den Komplementarstrahl ist 
also die ran Ac; enflang der Rotationsachse verschobene, durch die Winkelposition X x 
auf der Helix 17 definierte Strahlenquellenposition. 
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Der Komplementarstrahl, der von der um Aq verschobenen Strahlenquellenposition des 
Naherungsstrahles ausgeht, muss die Zylinderoberflache an der z-Position treffen, von 
der der dazugehorige direkte Strahl ausgeht Dies ergibt folgende Bestimmungs- 
gleichung fiir die Koordinate v x des Komplementarstrahles : 



tt.-tt.^.assL. ( 9) 



Daraus folgt durch Einsetzen von Gleichung (3) 



v, =<: 



*(2P 



/»(2p 



2i?cosp 

+ 

2 J? cos P 



fiir v 0 >0 



fiir v 0 <0 



(10) 



10 



Daher ergibt sich der Komplementarstrahl aus dem Naherungs strahl durch Verschieben 
der Strahlenquelle um Ac, entlang der Rotationsachse und durch Verschieben des 
Auftreffortes des Strahles auf dem Detektor um 



15 Av = v 1 -v 1 =-Aq— — L. (11) 

Die u-Koordinate des Komplementarstrahles ist gleich der entsprechenden Koordinate 
des Naherungsstrahles. 

20 Nachdem der Verlauf des Komplementarstrahles fiir einen direkten Strahl bekannt ist, 
kann der dem Komplementarstrahl zugeordnete Komplementarmesswert im Schritt 109 
bestimmt werden. Wenn der Komplementarmesswert akquiriert worden ist, kann der 
akquirierte Messwert als Komplementarmesswert verwendet werden. Dies ist aber in 
der Regel nicht der Fall, da meistens die Strahlenquellenposition des Komplementar- 

25 strahles nicht auf der Helix 17 liegt 
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10 



Ans „Improved 2D rebinning of helical cone-beam CT data using John's equation", 
Defrise, Noo, Kudo, Ml 0-74 in 2002 IEEE Nuclear Science Symposium Conference 
Record, guest editor Scott Metzler, 10-16 November, Norfolk, Virginia, USA, ISBN 0- 
7803-7637-4 ist bekannt, dass mit BBlfe von Johns-Gleichung ein zu den gemessenen 
Werten in guter Naherung konsistenter Komplementaimesswert durch Addition eines 
Korrekturmesswertes auf den Naherungmesswert berechnet werden kann. Der 
Korrektunnesswert AffC^u^vi) ist durch folgende Gleichung gegeben: 

^■ y -"fH^ t aiV l- t '" ) ^r ) l • ™ 



Zur Ermittlung des Korrekturmesswertes anhand Gleichung (12) konnen erfindungs- 
gemaB alle oder nur ein Teil der akquirierten Messwerte herangezogen werden. In 
diesem Ausfuhrungsbeispiel werden ausschliefilich Messwerte verwendet, deren 

15 zugeordnete Strahlen von Strahlenquellenpositionen y(X) ausgehen, die auf der Helix 
17 in der Nahe der Strahlenquellenposition y(X x ) des Naherungsstrahls angeordnet 
sind. Das heifit, es werden nur Messwerte verwendet, deren zugeordnete Strahlen- 
quellenpositionen im Bereich y(A,) mit X t - A x < X < X x + A x liegen. Hierbei ist A x so 
gewahlt, dass zehn unterschiedliche Strahlenquellenpositionen y(X) beriicksichtigt 

20 werden- In anderen Ausfuhrungsfonnen kann eine andere Anzahl an Strahlenquellen- 
positionen und damit an Messwerten verwendet werden. 

Der Bereich A x wird mindestens so groB gewahlt, dass redundante Messwerte zumin- 
dest zu einigen Komplementarmesswerten beitragen. Diese Bedingung ist allerdings in 
25 der Regel schon dadurch erfullt, dass der Naherungsmesswert selber redundant ist 

Zur Ermittlung des Korrekturmesswertes werden gemaB Gleichung (12) zunachst alle 
Messwerte, deren zugeordnete Strahlen von den Strahlenquellenpositionen ausgehen, 
die in dem oben genannten Bereich liegen, zum Erzeugen eines ersten Ableitungsdaten- 
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satzes partiell nach der Winkelposition X und der v-Koordinate abgeleitet Diese und 
auch die folgenden Ableitungen konnen bspw. mit der Methode der finiten Differenzen 
durchgefuhrt werden. AuBerdem werden die uispriinglichen, also die nicht abgeleiteten 
Messwerte, deren zugeordnete Strahlen von Sttahlenquellenpositionen ausgehen, die in 
5 dem oben genannten Bereich liegen, zum Erzeugen eines zweiten Ableitungsdaten- 
satzes zweimal partiell nach v abgeleitet. 

Die Messwerte des zweiten Ableitungsdatensatzes werden dann mit einem Wichtungs- 
faktor multipliziert, der quadratisch in u und linear in v ist In diesem Ausfiihrungs- 
1 0 beispiel ist dieser Wichtungsfaktor gleich 



Als nachstes werden die Messwerte des ersten und zweiten Ableitungsdatensatzes zu 
1 5 einem Satz von Messwerten addiert, und diejenigen resultierenden Werte, deren 
zugeordnete Strahlen von der Winkelposition X x ausgehen, auf die Detektorzeile vj 
und auf Spalten treffen, deren Koordinate u kleiner gleich u x ist, werden aufaddiert 

SchlieBlich wird die resultierende Summe mit einem Faktor multipliziert, der propor- 
20 tional zu Ac; und umgekehrt proportional zu R ist Bevorzugt wird die Summe mit 

Aq/R multipliziert. Die mit dem Faktor multiplizierte Summe ist der Korrekturmess- 
wert fur den Naherungsmesswert g(X t ,u l9 v x ) bzw. fur den Messwert g(X 0 , u 0 , v 0 ) des 
direkten Strahls. 

25 Nach Berechnung des Korrekturmesswert kaim der Komplementarmesswert 
g*(K> u Q> v *) zum Messwert g0 0 >*o> v o) 



Ruv-h(R 2 +u 2 ) 
R 2 



(13) 



g c (K> u o> v o) = g&>l>Ul>Vl) + 4ff(ki,«i, Vj) 



(14) 



30 durch Addition des Korrekturmesswertes auf den Naherungsmesswert bestimmt werden. 
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Nachdem die Schritte 105 bis 109 fur alle Messwerte g(X 0 , u 0 ,v 0 ) durchgefiihrt 
worden sind, die im PI-Fenster liegen, sind zwei Satze von Messwerten vorhanden, der 
im Schritt 103 akquirierte und der dazu komplementare Datensatz. 

5 

Im Schritt 111 ward jeder Messwert, fur den ein Komplementarmesswert ermittelt 
worden ist, durch eine Summe auf gewichtetem Komplementarmesswert und dazu- 
gehorigem Messwert ersetzt. Dazu wird der jeweilige Messwert auf den dazugehorigen 
Komplementarmesswert addiert, wobei vor der Addition jeder Wert mit einem Wich- 

10 tungsfaktor multipliziert wird Ist der Komplementarmesswert akquiriert worden, so 
werden dieser Komplementarmesswert und der dazugehorige Messwert gleich 
gewichtet^ insbesondere jeweils mit 0.5 multipliziert. Ist der Komplementarmesswert 
mit EQlfe einer Konsistenzbedingung ermittelt worden, bspw. mit Johns-Gleichung, so 
wird der Komplementarmesswert mit einem kleineren Wichtungsfaktor multipliziert als 

15 der dazugehorige Messwert. So kann der Messwert z.B. mit 0.9 oder 0. 8 multipliziert 
werden, wahrend der Komplementarmesswert mit 0. 1 oder 0.2 multipliziert wird. In 
diesem Ausfuhrungsbeispiel ist die Summe der Wichtungsfaktoren fur ein Messwert- 
Komplementarmesswert-Paar gleich 1. 

20 Im folgenden Schritt 1 13 wird ein CT-Bild, in diesem Ausfuhrungsbeispiel die 

Absorptionsverteilung im Untersuchungsbereich, aus den in den vorherigen Schritten 
ermittelten Messwerten, die im PI-Fenster liegen, also insbesondere aus den durch 
Addition von gewichteten Komplementarmesswerten und dazugehorigen Messwerten 
erzeugten Messwerten, mit einem exakten Verfahren rekonstruiert. Die einzelnen 

25 Rekonstruktionsschritte sind in Fig. 7 dargestellt 

In anderen Ausfiihrungsformen konnte auch nur fur einen Teil der Messwerte, die in 
einem Rekonstruktionsfenster liegen, ein Komplementarmesswert ermittelt worden 
sein. Dann wurden zur Rekonstruktion zum einen die im Rekonstruktionsfenster 
30 liegenden, durch eine Addition von gewichtetem Komplementarmesswert und 
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gewichtetem Messwert erzeugten Messwerte und zum anderen, falls fur einige im 
Rekonstruktionsfenster liegenden Messwerte kein Komplementarmesswert ermittelt 
worden ist, der jeweilige akquirierte Messwert verwendet werden. 

5 Zum Verstandnis der exakten Rekonstruktion wird zunachst folgende Gleichung aus El 
zitiert: 

2tc 2 J x) \x-yQ$isui! dq '*= x 

10 Diese Gleichung beschreibt eine exakte Rekonstruktion der Absorption durch Riick- 
projektion der Messwerte. Hier bezeichnet /(x) die raumliche Absorptionsverteilung 
im Untersuchungsbereich an der Stelle x. 

Der Messwert D f (y,&) lasst sich durch folgendes Linienintegral beschreiben: 
15 D f (y(X),®)= Jd//(y+/0) (16) 

0 

Der Einheitsvektor 0 gibt hierbei die Richtung des zum Messwert gehorenden Strahles 
an. 

Im Schritt 201 werden die Messwerte gemaB Gleichung (15) partiell nach q, d.h. nach 
20 der Winkelposition der Strahlenquelle an der Stelle q = X abgeleitet. Dabei ist zu 

beachten, dass nur y von q abhangt und nicht© , so dass fur die Ableitung jeweils Mess- 
werte von parallelen Strahlen beriicksichtigt werden miissen. Da paraDele Strahlen den 
gleichen Kegelwinkel aufweisen, treffen bei dem hier verwendeten fokus-zentrierten 
Detektor 16, wie in Fig. 8 dargestellt, die parallelen Strahlen 51 auf die gleiche Detek- 
25 torzeile 53. Hierbei ist der Kegelwinkel eines Strahles der Winkel, den dieser Strahl mit 
einer zur Rotationsachse 14 senkrechten Ebene einschlieBt. Zum partiellen Ableiten 
konnen die Messwerte zunachst umsortiert werden. Dazu werden Messwerte, die zu 
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parallelen Strahlen gehoren, also zu derselben Detektorzeile 53 aber zu unterschied- 
lichenWinkelpositionen X a , X b , X c der Strahlenquelle, jeweils zu einer Menge 
zusammengefesst Die Messwerte jeder Menge werden z. B. numerisch mit der 
bekannten Methode der finiten Differenzen nach der Winkelposition der Strahlenquelle 
5 abgeleitet 

Der Einheitsvektor © hangt vom K-Winkel y ab, der mit Hilfe sogenannter K-Ebenen 
55 beschrieben werden kann (siehe Fig. 9). Die K-Ebenen 55 werden im Folgenden 
erlautert 



10 



Zur Bestimmung einer K-Ebene 55 wird eine Funktion 



n (17) 

, A. > A>2 > A, — 2lt 



eingefuhrt, die von nicht-negativen, ganzzahligen Werten n und m, n>m, abhangt. In 
15 diesem Ausfuhrungsbeispiel wird n - 2 und m=l gewahlt Es konnten aber auch 
andere Werte n, m gewahlt werden. Die Gleichung (15) wiirde trotzdem exakt bleiben, 
nur die Lage der K-Ebenen 55 wiirde sich andern. Des weiteren werden die 
Vektorfunktion 



20 



u(X^X 2 ) = 



|y(^(^A 2 ))-y(X)Jx|v(X 2 )-y(X)] Sgn(X2 - X) ' 0< k" 5 l 



y(X)xy(X) 



(18) 



Jy(^)xy(M 3 
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und der Einheitsvektor 



P(^,x) = 



|x-y(M 



(19) 



5 definiert Der Vektor p zeigt von der Strahlenquellenposition y(X) zur Position x. Zur 
Bestimmung der K-Ebene wird nun ein Wert X 2 e I n (x) gewahlt, so dass yX, 
ytynXM)), y(^2) und x in einer Ebene liegen. Diese Ebene wird als K-Ebene 55 und die 
Schnittlinie zwiscben der K-Ebene 55 und der Detektorflache 74 wird als K-Linie 57 
bezeichnet. Der Detektor 74 ist in Fig. 9 durch zwei aufeinander folgende Windungen 

10 der helixrormigen Trajektorie 1 7 begrenzt und weist die Kriimmung der Helix 1 7 auf. 
Dieser Detektor wurde hier als Beispiel zur Veranschaulichung genutzt. Fiir andere 
Detektoren, wie den fokus-zentrierten oder den planaren Detektor, konnen entsprechen- 
de Schnittlinien 57 bestimmt werden. In Fig. 9 ist ein facherformiger Teil einer K-Ebene 
dargestellt. Die Rander des Fachers treflen sich am Strahlenquellenort. Diese Definition 

15 der K-Ebene 55 ist aquivalent zur Losung der Gleichung 



nach %2- Damit ist u der Normalenvektor der K-Ebene 55. Zur Bestimmung der 
20 Vektorfunktion® (X, x, y) wird der Vektor 



definiert Mit den Definitionen fur p und e kann nun die Vektorfunktion ®(X, x, y) wie 
25 folgt angegeben werden: 



(x - y(X)) -u(X,X 2 ) = 0, X 2 e /, 



(20) 



e(X,x) = P(X,x)xu(X,x) 



(21) 



®(X, x, y) = cos y • Q(X, x) + sin y • e(X, x) 



(22) 
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Da beide Vektoren, 0 und e, senkrecht zu u orientiert sind, gibt der K-Winkel y die 

Richtung des Vektors Qjund damit die Richtung eines Strahls, innerhalb einer K-Ebene 
an. 

5 

Die K-Ebenen und K-Linien sind im Einzelnen in El beschrieben, worauf hiermit Bezug 
genommen wird. 

Im Schritt 203 werden gemaB Gleichung (15) die abgeleiteten Messwerte entlang von k- 
10 Linien mit einem Wichtungsfaktor multipliziert, der mit zunehmendem Sinus des k- 
Winkels yabnimmt und insbesondere gleich dem Kehrwert des Sinus des K-Winkels ist, 
und aufaddiert. Dazu wird fiir jede Stelle x im Untersuchungsbereich und fur jede 
Strahlenquellenposition X eine K-Linie bestimmt, wobei, wie oben erlautert, ein Wert 
X 2 e 7 PI (x) so ausgewahlt wird, dass y(X), y(kAA, 2 )), x ^ einer Ebene liegen, 

15 der K-Ebene. Die K-Linie wird dann als Schnittlinie zwischen der K-Ebene und der, 
Detektorflache ermittelt Die Multiplikationen mit dem Wichtungsfektor und die 
Integrationen bzw. Additionen konnen beispielsweise mit EQlfe einer Fouriertrans- 
formation durchgefuhrt werden. 

20 Die abgeleiteten und integrierten Messwerte lassen sich durch folgende Gleichung 
darstellen: 

p(y(X),*(X,x» = )JLJL Df (j( q ),®(X,x>Y)) U ( 23 ) 
_;siny oq H 

Hier bezeichnen p(y(X), <I> (X, x)) die abgeleiteten und integrierten Messwerte und 
25 ® (X, x) einen Einheitsvektor, der von der Strahlenquellenposition y(X) in Richtung der 
Stelle x im Untersuchungsbereich zeigt. 
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Der noch fehlende Integrationsschritt in Gleichung (15) bzw. die Ruckprojektion der 
Messwerte kann nun durch folgende Gleichung beschrieben werden: 

5 

Vor der Ruckprojektion kann im Schritt 207 ein Rebinning der Messwerte erfolgen. 
Durch das Rebinning werden die Messwerte so umsortiert und uminterpoliert, als waxen 
sie mit einer anderen Sttahlenquelle (einer ausgedehnten, auf einem Teil einer Helix 
angeordneten Strahlenquelle, die jeweils zueinander parallele Strahlenfacher emittieren 
1 0 kann) gemessen worden. 

Dies wird anhand von Fig. 10 naher erlautert Mit 17 ist dabei die helixfbrmige Trajek- 
torie bezeichnet, von der aus die Strahlenquelle den Untersuchungsbereich durchstrahlt. 
Mit 43 ist ein facherfonniges Strahlenbundel bezeichnet, dass von der Strahlenquellen- 
position S 0 ausgeht und dessen Strahlen in einer die Rotationsachse 14 enthaltenden 
15 Ebene verlaufen. Man kann sich das kegelformige Strahlenbundel, das von der Strahlen- 
quelle in der Position S 0 emittiert wird, aus einer Vielzahl von ebenen Strahlenfachern 
zusammengesetzt denken, die sich in zur Rotationsachse 14 parallelen Ebenen befinden 
und sich in der Strahlenquellenposition S Q schneiden. Fig. 10 zeigt von diesen Strahlen- 
fachern nur einen einzigen, namlich den Strahlenfacher 43. 

20 

AuBerdem sind in Fig. 10 noch weitere Strahlenfacher 41, 42 und 44, 45 dargestellt, die 
parallel zu dem Strahlenfacher 43 sind und in zueinander und zur Rotationsachse 14 
parallelen Ebenen liegen. Die zugehorigen Strahlenquellenpositionen S_ 2 , S_{ und 5, , 
werden von der Strahlenquelle S eingenommen, bevor bzw. nachdem sie die 
25 Strahlenquellenposition S 0 erreicht hat Alle Strahlen in den Strahlenfachern 41 bis 45 
haben denselben Projektionswinkel. Die StrahlenScher 41 bis 45 definieren ein 
Strahlenbiindel 70 mit einer zeltartigen Form. 
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Die nach dem Rebinniag ermittelten Messwerte werden anschliefiend zurRekonstruk- 
tion der Absoiptionsverteilung im Untersuchungsbereich dutch eine Ruckprojektion, in 
diesem Ausfuhrungsbeispiel gemaB Gleichung (23), herangezogen. 

5 Dazu wird im Schritt 209 ein Voxel F(x) innerhalb eines vorgebbaren Bereiches (field 
of view — FOV) im Untersuchungsbereich und ein Pl-Ihtervall I n (x) fur diesen Voxel 
bestimmt. Dann wird im Schritt 21 1 eine Winkelposition X innerhalb des Ihtervalls 
^pi ( x ) vorgegeben. Im Schritt 213 wird gepriift, ob fur die Winkelposition X ein 
Messwert vorhanden ist, dessen Strahl durch die Mitte des Voxels V(±) verlauft Kaon 

10 ein solcher Strahl nicht gefunden werden, so wird ermittelt, an welcher Stelle ein 

mittdger Strahl auf die Detektorflache getroffen ware. Der dazugehorige Messwert. wird 
dann durch Interpolation der Messwerte von benachbarten Strahlen berechnet. Der 
Messwert, der dem den Voxel passierenden Strahl zugeordnet werden kann, bzw. der 
durch Interpolation gewonnene Messwert wird im Schritt 214 mit einem Wichtungs- 

15 faktor multipliziert, der mit zunehmenden Abstand der Strahlenquelle y(X) von dem zu 
rekonstruierenden Ort x im Untersuchungsbereich Heiner wird In dieser Ausfuhrungs- 
fonn ist dieser Wichtungsfektor gemaB Gleichung (24) gleich 1 /|x - y (X)\ . Im Schritt 
215 wird der gewichtete Messwert auf den Voxel V(x) akkumuliert. Im Schritt 217 
wird gepriift, ob alle Winkelpositionen X im Intervall 7 PI (x) betrachtet worden sind. Ist 

20 dies nicht der Fall, so verzweigt das Ablaufdiagramm zum Schritt 21 1 . Ansonsten wird 
im Schritt 219 gepriift, ob alle Voxel V(±) im FOV durchlaufen sind. 1st dies nicht der 
Fall, so wird mit Schritt 209 fortgefahren. Wenn dagegen alle Voxel V(±) im FOV 
durchlaufen sind, so ist die Absorption im gesamten FOV und damit das CT-Bild 
ermittelt worden, und die exakte Rekonstruktion des CT-Bildes ist im Schritt 221 

25 beendet 

Wie oben bereits erwahnt, kann sich die Strahlenquelle in anderen Ausfuhrungsformen 
relativ zum Untersuchungsbereich auch auf einer kreisformigen Trajektorie bewegen. 
Wenn in einem solchen Fall zur Rekonstruktion eines CT-Bildes nur Messwerte 
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verwendet werden, die in einem Rekonstruktionsfenster liegen, so konnen auch hier 
zumindest fur einige Messwerte Komplementarmesswerte insbesondere mit Johns- 
Gleichung eimittelt werden. Aus einem Messwert-Komplementarmesswert-Paar kann, 
wie oben fur die helixfonnige Trajektorie beschrieben, ein Summe gebildet werden, die 
5 den Messwert im akquirierten Datensatz ersetzt Aus den so erzeugten Daten kann 
schlieBlich mit bekannten Rekonstruldionsverfahren ein CT-Bild rekonstruiert wird. 

Ein Computertomographie-Verfahren, bei dem sich die Strahlenquelle relativ zum 
Untersuchungsbereich auf einer kreisfiiixnigen Trajektorie bewegt und ausschliefllich 
10 Messwerte verwendet werden, die in einem Rekonstruktionsfenster liegen, ist bspw. aus 
„A fast and efficient method for sequential cone-beam CT" Koehler, Proksa, Grass, 
Medical Physics, vol. 28, no. 11, pp. 2318-2327 bekannt. 
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BEZUGSZEICHENLISTE 

OfEnungswinkel 

X, X, , X 2 , % a ,X h , X 0 Winkelpositionen der Strahlenquelle auf der helixformigen 



> 




irajeKtone 






Verschiebung in z-Richtung 




Av 


Verschiebung in v-Richtung 






Stelle im Untersuchungsbereicb 






Helixabschnitt 


10 


s 


Strahlenquelle 






Strahlenquellenpositionen 




1 


Gantry 




2,5 


Motor 


15 


3 


Kollimatoranordniin or 

Jfcg^yr Hill M Hp1>%#JL. MAJ 11 "It 1 




4 


Strahlenbiindel 




7 


Steuereinheit 




10 


Bildverarbeitunesrechner 




11 


Monitor 


20 


13 


Untersuchungsbereich 




14 


Rotationsachse 




16 


Detektoreinheit 




17 


helixfiirmige Trajektorie 




21,23 


PI-Linien 


25 


25 


PI-Fenster 




31 


Naherungsstrahl 




33 


direkter Strahl 




35 


Zylinderoberflache 




37 


PI-Gerade 


30 


41, 42, 43, 44, 45 


zueinander parallele Strahlenfacher 




51 


parallele Strahlen 




53 


Detektorzeile 




55 


K-Ebene 




57 


K-Linie 


35 


60 


planarer Detektor 




70 


Strahlenbiindel 




74 


Detektorflache 
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PATENTANSPRftCHE 



1 . Computertomographie-Verfehren mit den Schritten: 

a) Erzeugen eines kegelformigen, einen Untersuchungsbereich (1 3) und ein darin 
befindliches Objekt durchsetzenden Strahlenbiindels (4) mit einer Strahlenquelle 
(S), 

b) Erzeugen einer Relativbewegung zwischen der Strahlenquelle (S) einerseits und 
dem Untersuchungsbereich (13) andererseits, die mindestens eine Rotations- 
bewegung um eine Rotationsachse (14) umfasst und insbesondere die Form 
einer Helix (1 7) oder eines Kreises hat, 

c) Akqumeren von Messwert^ die von der 

jenseits des Untersuchungsbereichs (13) abhangen, mit einer Detektoreinheit 
(1 6) wahrend der Relativbewegung, 

d) Ermitteln von jeweils einem Komplementarmesswert fur zumindest einige der in 
Schritt c) akquirierten und in einem Rekonstruktionsfenster (25) liegenden 
Messwerte mit Hilfe von redundanten Messwerten, wobei die dem jeweiligen 
Messwert und dem dazugehorigen Komplementarmesswert zugeordneten 
Strahlen entgegengesetzt zueinander orientiert sind, 

e) Ersetzen jedes Messwertes, fur den in Schritt d) ein Komplementarmesswert 
ermittelt worden ist, durch eine Summe aus gewichtetem Messwert und 
gewichtetem Komplementarmesswert, 

f) Rekonstruieren eines CT-Bildes des Untersuchungsbereichs (1 3) aus den im 
Rekonstruktionsfenster (25) liegenden Messwerten. 
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2. Computertomograpme-Verfahren nach Anspruch 1, 
dadurch peke rmgftihhnet, 

5 dass, wenn ein Strahl, der einem Komplementarmesswert nach Schritt d) zugeordnet ist, 
gleich einem Strahl, der einem im Schritt c) akquirierten Messwert zugeordnet ist, 
verlauft, der Komplementarmesswert gleich diesem Messwert gesetzt wird, und dass, 
wenn ein Strahl, der einem Komplementarmesswert nach Schritt d) zugeordnet ist, nicht 
gleich einem der Strahlen, die den im Schritt c) akquirierten Messwerten zugeordnet 

1 0 sind, verlauft, dieser Komplementarmesswert mit Hilfe von Johns-Gleichung ermittelt 
wird. 

3. Computertomographie-Verfahren nach Anspruch 1, 
dadurch gekennzeichnet 

1 5 dass vor der Addition im Schritt e) der Komplementarmesswert und der dazugehorige 
Messwert mit jeweils einem Wichtungsfaktor multipliziert werden, wohei die 
Wichtungsfaktoren fur einen Komplementarmesswert und einen dazugehorigen 
Messwert gleich sind, wenn der Strahl, der dem Komplementarmesswert nach Schritt d) 
zugeordnet ist, gleich einem der Strahlen, die den im Schritt c) akquirierten Messwerten 

20 zugeordnet sind, verlauft, und der Wichtungsfaktor fur einen Messwert grofier ist als der 
Wichtungsfaktor eines dazugehSrigen Komplementarmesswerts, wenn der Strahl, der 
dem Komplementarmesswert nach Schritt d) zugeordnet ist, nicht gleich einem der 
Strahlen, die den im Schritt c) akquirierten Messwerten zugeordnet sind, verlauft. 



25 
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4. Computertomographie-Verfahren nach Anspruch 1, 
dadurch ft ekenn zeichnet 

dass im Schritt b) die Relativbewegung die Form einer Helix (1 7) hat und dass das 
Rekonstruieren eines CT-Bildes im Schritt f) folgende Schritte aufweist: 

- partielles Ableiten von Messwerten, denen parallele Strahlen (51) mit 
unterschiedlichen StrahlenqueUienpositionen ( X n ,X h , X e ) zugeordnet sind, nach 

einer dem jeweiligen Messwert zugeordneten Winkelposition der Strahlenquelle 
(S)aufder Helix (17), 

- Filtem der abgeleiteten Messwerte entlang von K-Linien (57), 

- Rekonstruktion des CT-Bildes durch Ruckprojektion derjenigen gefilterten 
Messwerte, die in einem PI-Fenster (25) liegen. 

5. Computertomographie-Verfahren nach Anspruch 4, 

dadurch ff ekennrairJitiftt 

dass das Filtem eines Messwertes folgende Schritte auj^veist: 

- Bestimmen einer K-Linie (57) fur den Messwert, 

- Multiplizieren derjenigen Messwerte, die auf der K-Linie (57) liegen mit jeweils 
einem Wichtungsfaktor, der mit dem Kehrwert des Sinus des k- Winkels 
zunimmt, insbesondere gleich diesem Kehrwert ist, 

- Aufeddieren der auf der K-Linie (57) liegenden, gewichteten Messwerte, wobei 
die resultierende Summe der gefilterte Messwert ist. 

6. Computertomograph mit 

- einer Strahlenquelle (S) zum Erzeugen eines kegelfiJrmigen einen 
Unteisuchungsbereich (13) und ein darin befindliches Objekt durchsetzenden 
Strahlenbiindels (41), 

- einer Antriebsanordnung (2,5), urn ein im Untersuchungsbereich (1 3) 
enthaltenes Objekt und die Strahlenquelle (S) relativ zueinander urn eine 
Rotationsachse (14) rotieren und sich parallel zur Rotationsachse (14) 
verschieben lassen zu konnen, 
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einer mit der Strahlenquelle (S) gekoppelten Detelctoreinheit (16), die eine 

Detektorflache aufweist, zur Akquisition von Messwerten, 

einer Rekonstniktionseinheit (10) zur Rekonstruktion des CT-Bildes innerhalb 

des Untersuchungsbereichs aus den von der Detektoreinheit (16) akquirierten 

Messwerten, 

einer Steuereinheit (7) zur Steuerung der Strahlenquelle (S), der Detektoreinheit 
(16), der Antriebsanordung (2,5) und der Rekonstruktionseinheit (10) 
entsprechend den folgenden Schritten: 

a) Erzeugen eines kegelformigen, den Untersuchungshereich (13) und das darin 
befindliches Objekt durchsetzenden Strahlenbundels (4) mit der 
Strahlenquelle (S), 

b) Erzeugen einer Relativbewegung zwischen der Strahlenquelle (S) einerseits 
und dem Untersuchungshereich (13) andererseits, die mindestens eine 
Rotationsbewegung urn die Rotations achse (14) umfasst und insbesondere 
die Form einer Helix (17) oder eines Kreises hat, 

c) Akquirieren von Messwerten, die von der Intensitat in dem Strahlenbundel 
(4) jenseits des Untersuchungsbereichs (13) abhangen, mit der 
Detektoreinheit (16) wahrend der Relativbewegung, 

d) Ermitteln von jeweils einem Komplementarmes swert fur zumindest einige 
der in Schritt c) akquirierten und in einem Rekonstruktionsfenster (25) 
liegenden Messwerte mit Hilfe von redundanten Messwerten, wobei die dem 
jeweiligen Messwert und dem dazugehorigen Komplementarmesswert 
zugeordneten Strahlen entgegengesetzt zueinander orientiert sind, 

e) Ersetzen jedes Messwertes, fur den in Schritt d) ein Komplementarmesswert 
ermittelt worden ist, durch eine Summe aus gewichtetem Messwert und 
gewichtetem Komplementarmesswert, 

f) Rekonstruieren eines CT-Bildes des Untersuchungsbereichs (1 3) aus den im 
Rekonstruktionsfenster (25) liegenden Messwerten. 
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7. Computerprogramm fiir eine Steuereinheit (7) zur Steuerung einer Strahlenquelle (S) 3 
einer Detektoreinheit (16), einer Antriebsanordnung (2,5) und einer Rekonstruktions- 
einheit (10) eines Computertomographen zur Durchffihrung des Verfahrens nach 
Anspruch 1 gemaB folgendem Ablauf: 

a) Erzeugen eines kegelffirmigen, einen Untersuchungsbereich (13) und ein 
darin befindliches Objekt durchsetzenden Strahlenbundels (4) mit der 
Strahlenquelle (S), 

b) Erzeugen einer Relativbewegung zwischen der Strahlenquelle (S) einerseits 
und dem Untersuchungsbereich (13) andererseits, die mindestens eine 
Rotationsbewegung urn eine Rotationsachse (14) umfosst und insbesondere 
die Form einer Helix (17) oder eines Kreises hat, 

c) AkqumerenvonMesswerten,die™^ 

(4) jenseits des Untersuchungsbereichs (13) abhangen, mit der 
Detektoreinheit (1 6) wahrend der Relativbewegung, 

d) Ermitteln von jeweils einem Komplementaimesswert fur zumindest einige 
der in Schritt c) akquirierten und in einem Rekonstruktionsfenster (25) 
liegenden Messwerte mit Hilfe von redundanten Messwerten, wobei die dem 
jeweiligen Messwert und dem dazugehorigen Komplementarmesswert 
zugeordneten Strahlen entgegengesetzt zueinander orientiert sind, 

e) Ersetzen jedes Messwertes, fur den in Schritt d) ein Komplementarmesswert 
ermittelt worden ist, durch eine Summe aus gewichtetem Messwert und 
gewichtetem Komplementarmesswert, 

f) Rekonstruieren eines CT-Bildes des Untersuchungsbereichs (13) aus den im 
Rekonstruktionsfenster (25) liegenden Messwerten. 
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ZUSAMMENF AS SUNG 

ComputertomograpMe-Verfahren mit redundanten Messwerten 

Die Erfindung betrifft ein Computertomographie-Verfehren, bei dem sich eine Strahlen- 
quelle relativ zu einem Untersuchungsbereich insbesondere auf einer helixfonnigen 
oder kreisfbrmigen Trajektorie bewegt. Messwerte werden mit einer Detektoreinheit 
akquiriert, wobei aus diesen Messwerten ein CT-Bild des Untersuchungsbereichs 
rekonstruiert wird. Bei der Rekonstruktion wird zumindest fur einige Messwerte, die in 
einem Rekonstruktionsfenster liegen, jeweils ein Komplementarmesswert bestimmt, 
dessen Strahl entgegengesetzt parallel zu dem Strahl des jeweiligen akquirierten 
Messwerts orientiert ist. Zur Berechnung der Komplementannesswerte, insbesondere 
mit EQlfe von Johns-Gleichung, werden redundante Messwerte verwendet Die Mess- 
werte, fur die Komplementannesswerte bestimmt worden sind, werden jeweils durch 
eine Summe aus gewichtetem Messwert und gewichtetem Komplementarmesswert 
ersetzt, und ein CT-Bild wird aus den ersetzten und gegebenenfalls aus akquirierten 
Messwerten, die im Rekonstruktionsfenster liegen, insbesondere mit einem exakten 
Rekonstruktionsverfaliren, rekonstruiert. 

Fig. 1 
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Fig. 3 




Fig. 4 
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